Medium modification of two-particle scattering in nonideal Bose systems 
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I. INTRODUCTION 

Since the recent observation of Bose-Einstein conden- 
sation (BEC) in rubidium 87 Rb 0, lithium 7 Li § and 
sodium 23 Na j| in atom traps at extremly low temper- 
atures new interest has arisen in the investigations of 
interacting Bose gases. For a comprehensive overview 
about the current experimental and theoretical activities 
see Q. The effect of finite size on the critical temper- 
ature of an ideal Bose gas, which is confined in a small 
cube, was investigated in However, when perform- 
ing experiments at higher densities, the question of the 
treatment of interaction in a nonideal Bose gas becomes 
essential. 

In |^] the temperature dependence of the moment of 
inertia of a dilute Bose gas in a harmonic trap and the 
role of interaction is discussed. A common attempt to 
include interaction is the approximation of a hard sphere 
Bose gas which can be treated analytically (To). There, 
the interaction can be replaced by a pseudopotential con- 
taining the free scattering length ao as the only param- 
eter. For the atoms relevant in the recent BEC exper- 
iments, values for the scattering length are obtained in 
the low density limit from measurements of the cross sec- 
tion jllj or from evaluation of spectroscopic data p2|,^3| . 
However, the isolated two-particle scattering properties 
should be modified in dense systems due to the medium 
as shown below. This means that the scattering processes 
in a dense medium cannot be interpreted in terms of the 
free scattering length only. 

The medium modifications of two-particle properties 
result in different macroscopic effects. First, nonequilib- 
rium processes are influenced. Usually, kinetic equations 
are derived in RPA or Born approximation jf^-Q]. The 
corresponding collision terms include final state Bose en- 
hancement. It was demonstrated [Q, that approaching 
the critical density of Bose condensation the scattering 
rates are increased due to Bose enhancement. In addi- 
tion, the influence of the medium can be included by a 
meanfield approximation in the framework of the quasi- 



particle concept. Employing the Bogoliubov transfor- 
mation, a quasiparticle approach was proposed by H,[|, 
which is applicable also below the critical temperature. 

However, additional effects arise if we go beyond 
a quasiparticle Born approximation as shown below. 
These effects are obtained by the solution of the Bcthe- 
Salpeter equation containing intermediate Bose enhance- 
ment. They cause a modification of the T-matrix or 
the corresponding transition matrix elements which are 
absent in Born approximation. This T-matrix leads to 
medium dependend effective scattering phase shifts and 
cross sections which are an important ingredient for col- 
lision integrals. In this paper we will not discuss kinetic 
theory in more detail (see p"§|]). 

Second, the thermodynamic quantities are modified 
by medium dependend two-particle properties. As well- 
known, nonideality effects of thermodynamic quantities 
can be expressed in terms of two-particle scattering phase 
shifts. A relation between these phase shifts and the sec- 
ond virial coefficient can be established using the Beth- 
Uhlenbeck formula iQ. It is possible to generalize the 
treatment of the isolated two-particle problem to dense 
quantum gases. This is achieved using the T-matrix with 
intermediate Bose enhancement as we will show below. 

Systematic attempts for the derivation of second virial 
corrections for dense quantum gases are given using ther- 
modynamic Green function technique for Fermi systems 
|f5| , ^6| or Ursell-operator technique [jl7| . In this paper we 
apply thermodynamic Green functions for Bose systems 
and study the microscopic two-particle scattering prob- 
lem in medium and the resulting equation of state. We 
discuss the conclusions for two examples with repulsive 
interaction, sodium 23 Na and rubidium 87 Rb, which play 
a major role in the BEC experiments. 

First, let us remind some definitions and results of the 
scattering of two interacting particles in vacuum. The 
free scattering length ao is defined in the limit of vanish- 
ing relative momentum k 
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ciq = — lim 



S (k) 



k^O k 
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This follows from the phase shift So(k) between the 
asymptotic incoming and outgoing plane wave of the 
scattered particle, which can be expanded for small rela- 
tive momentum 



cot<Wfe) = l -r + \kr Q +0(k 3 ) 

Clgk 2 



(2) 



The coefficient tq in the second term characterizes the 
range of the potential. 

Let us consider an imperfect Bose gas with <C A 
where A = ( 2 Z h r r ) 1 ^ 2 is the thermal wavelength. In 
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this case, a zero range pseudopotential can be employed 
which is parametrized by the free scattering length ao 
only. For densities below the critical density of the ideal 
Bose-Einstein condensation tibec = 2.612A~ 3 , a first or- 
der correction of the ideal Bose gas pressure pid can be 
derived [Hoi 



p(n, T) = Pid(n, T) H n 



(3) 



In the nondegenerate limit (n -C bbec), the equation of 
state (j^) can be rewritten in terms of a virial expansion 
either of pressure or of density 



p(n, T) — nksT + B p (T)n 2 kBT 
n((J,,T) = m<i((J,,T) + B n (T)n 2 d (fJ,,T) 



with the second virial coefficient B p (T) = 2ao\ 2 or 
B n {T) = — 4aoA 2 . In the following, we present a method 
of deriving the in-medium generalizations of the scatter- 
ing phase shifts and the scattering length as well as an 
equation of state for a dense Bose gas with respective 
virial corrections in the normal phase. 



II. GREEN FUNCTION APPROACH TO 
IN-MEDIUM SCATTERING AND 
THERMODYNAMICS 

In a dense medium, the single-particle, two-particle, 
etc., energies are modified compared to vacuum values. 
For instance, quasiparticle energies and exchange pro- 
cesses with the medium have to be included. A system- 
atic quantum statistical approach can be given employing 
thermodynamic Green functions. In the case of fermions, 
the in-medium scattering problem was investigated in 
||il|fl6l| . We follow the method as outlined in Jl6| and 
apply it for the case of bosons. 
We start from the Dyson equation 



G\ — Gi -\- G-^EGi 



(5) 



which relates the free one-particle Green function 
G5(l; z) — (z — ei)" 1 and the the full one-particle Green 
function G\ in a selfconsistent way. The latter one can 
be written as 



Gi(l;«) = (2-ei-E(l;«))- 1 



(6) 



where 1 = {p\ , o\ } denotes quantum numbers momen- 

h 2 k 2 

turn and spin and e\ = -^zr is the free single particle 
energy. 

The self energy S can be expressed in terms of the T- 
matrix and the spectral function A as 
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E ' 

2 

" + °° dE 

7T 



dE 2 
2tt 



A(2, E 2 



Im 7^(12, 12; E) 



51(^2)^(12,12) 
gi{E 2 )-g 2 {E) 



E-z-Ei 



(7) 



where the index "ex" denotes the sum of direct and ex- 
change matrix elements in the respective quantity and 
g n (E) = {exp((_E — n/i)/T) — l}^ 1 is the n-particle Bose 
distribution function. 

Using thermodynamic Green functions, the T-matrix 
of the quantum scattering theory is generalized to finite 
temperature and density. It contains a partial summation 
of the interaction. The resulting Bethe-Salpeter equation 
for the T-matrix within the ladder approximation reads 
111 



(4) T(12, 1'2'; z) = V(12, 1'2' 



V(12, 34) G^(34, 56; z) T(56, 1'2'; 



;uri<i 



(8) 



The T-matrix describes the two-particle properties and 
corresponds to the Schrodinger equation but contains the 
influence of the medium in the uncorrelated two-particle 
Green function G°. The uncorrelated two-boson Green 
function reads in spectral representation 

r+ °° dE 1 f +0 ° dE 2 l + g 1 (E 1 )+g 1 (E 2 ) 



2tt 



E\ — E 2 



xA(l,E 1 )A(2,E 2 )S n ,S 22 , 



(9) 



The one-particle spectral function is defined as A(l, E) — 
i[Gi(l,E + iQ) - Gx(l,E - iO)] and reads by using the 
Dyson equation (||) 



A(1,E) 



2ImE(l,J5) 



[E-ex -ReE(l,£;)] 2 + [ImS(l,£)] ; 



(10) 

The set of equations ([?]) , (@) , @ and (|o|) has to be solved 
selfconsistently. 
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The T-matrix is directly related to the two-particle 
scattering properties. E.g., in the continuum of scatter- 
ing states the phase shifts S can be expressed as 



cot <5(12, 12; E) 



ReT(12, 12; E) 
ImT(12,12;£) 



(11) 



which contain the influence of the medium. These in- 
medium scattering phase shifts are a direct generalization 
of the standard definition of free phase shifts which are 
included in the low density limit [ ^6[ . 

From the spectral function thermodynamic properties 
of the system can be derived. The one-particle density 
of a Bose system reads 



l ^' T ) = ^E / ^9i{E)A{l,E) 



(12) 



where f2 is a normalization volume. 



III. EXTENDED QUASI-PARTICLE 
APPROXIMATION AND GENERALIZED 
BETH-UHLENBECK FORMULA 

The full selfconsistent solution of eqs. (Q) to (|H 
is a complicated task. Therefore, we consider the 
spectral function (|l0|) in quasi-particle approximation 
A qp (l,E) = 2tt5{E - ei - A(l)) which shows delta- 
peaked spectral weight at sharp quasi-particle energies 
ei = ei + A(l) of the single bosons. The single-particle 
energy E\ is shifted by A(l) = Re S(l, z)\ z=ei which is 
the real part of the self energy at the quasi-particle 
energy. The quasi-particle spectral function is applied 
to the self energy (^h and to the uncorrelated two-boson 
Green function G° (a) which enters the ladder T-matrix 
(|J) . From this a spectral function can be obtained which 
goes beyond the quasi-particle approximation. The so- 
called extended quasi-particle approximation is obtained 
by assuming that the imaginary part of the self energy is 
small. The resulting spectral function can be represented 
as a 5-function and an additional contribution due to cor- 
related particles in form of a derivation of the principal 
value Il5l|il| 



A cqp (l,E) 



2nS(E-e 1 ) 



l-^Re£(M)U =ei 
d P 

2ImE(l,S)- 
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(13) 



Applying ([13]) in cq. (|l2|) results in an equation of state 
for the total density which contains quantum statisti- 
cal spectral weight from single bosons in quasi-free or in 
possible two-particle bound and scattering states. Simi- 
lar to the classical equation of Beth and Uhlenbeck , 
the total density splits into density contributions due to 
free quasi-particlcs and correlated particles in bound and 



scattering states. In an analogous way as done in [[L6| for 
a fermionic system we give a generalized Beth- Uhlenbeck 
formula for the total density in a bosonic system which 
is only valid in the normal phase for temperatures above 
the onset of Bose-Einstcin condensation 



«tot(/i, T) = nfree(M) T) + 2n corr (^, T) 



(14) 



with 



1 
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a { ^ K(K; E) - i sin(2<5 Q (K; E))] } 



x J2 C 



where, if necessary, summations over bound states and 
scattering channels a with a corresponding degeneracy 
factor c a have to be performed. The in-medium bound 
state energy E^(K) depends on total momentum K, 
chemical potential and temperature. Like the scatter- 
ing energy E, it is taken relative to the energy at the 
continuum edge E cont (K) = h 2 K 2 /Am + 2A(K) which 
separates bound state (E^(K) < 0) and scattering state 
region (E > 0). 
Rewriting ( |l4| ) as 

n to t(^,T) =n ilcc (p,T)[l + B*(n,T) n frcc (^,T)] (15) 

for (rtfrcc < tibec): we can define a quantity B*(fi,T) — 
2n. corr (/^, T)/(nf ree (A*, T)) 2 which represents the non- 
ideality contribution in analogy to the virial expansion of 
the total density (Q). We remark that B* contains non- 
ideality contributions due to two-particle interaction and 
degeneracy in any order of density, but no three-particle 
interaction contributions. While the second virial coeffi- 
cient B n (T) per definitionem contributes to second order 
in density and depends on temperature only, B* also de- 
pends on chemical potential and hence on density. In 
the low density and low temperature limit, or more ex- 
act for a <C A and rif rco A 3 <C 1, the non-ideality coeffi- 
cient B*(fi, T) coincides with the second virial coefficient 
B n (T) fromeq. (@). 



IV. MODEL CALCULATION 

In the case of a separable representation of the interac- 
tion, the Bethe-Salpeter equation (||) can be solved and 
the T-matrix can be written in an analytic form. The 
choice of a separable potential is motivated by the fact 
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that there is a rigorous method to replace an arbitrary 
potential by a rank- TV separable potential Sepa- 
rable potentials are applied in the description of differ- 
ent physical systems, e.g., by Nozieres and Schmitt-Rink 
p0| for electron-hole interaction or in nuclear physics for 
the nucleon-nucleon interaction |2l]] . As shown below the 
contact interaction which is widely used in atomic physics 
is a special case of a separable potential. 

For our model calculation we choose a rank-one sepa- 
rable potential of the form 



7(12, 1'2') = V(kK, k'K') = V v(k)v(k')S_ 



KK' 



with a formfactor 



v{k) 



k 2 + (3 2 



(16) 



(17) 



The potential depends on the square of the incoming and 
outgoing relative wave vector k = \ki — k 2 \/2 and k' = 
\k[ — k' 2 \/2 and the two parameters (3 (« inverse range) 
and Vq (strength), but it is independent of the center- 
of-mass momentum TiK. For j3 — > oo the special case 
of a contact potential, i.e. a zero range interaction, is 
included. 

Now, using relative coordinates, we can solve the in- 
medium Bethe-Salpeter equation (|^) explicitly 



T(kK, k'K';z) 



V{kK, k'K') 
1 - J{K;z) 



with 



J{K-z)= V(k'K,k'K')G° 2 (k'K; 



(18) 



(19) 



k'.K' 



The two-boson Green function in quasi-particle approxi- 
mation reads 



G° 2 (kK;z) = 



Q(K,k) 
' — ei - £2 



(20) 



The medium effects enter via the Bose enhancement fac- 
tor Q(K, k) and via the self energy shift of the quasi- 
particle energies. For our exploratory calculation below 
we replace the shift by its thermal average A. This rigid 
shift can be incorporated in an effective chemical poten- 
tial n* = fi — A. 

The integration of angles in the Bose distributions 
Q(K, k) = J f£ [1 + gi(ei) + .9i (e 2 )] gives an analytic ex- 
pression for the Bose enhancement factor 



Q(K, k) = 1 



2mfcs7 1 
h 2 Kk 



In 
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2mk B T 



k B T 



(21) 



Now, the quantity J reads explicitly 



J(K, E(k) + iO) = ^ I dk' k 



Vo_ f AU , , J2 v(k') 2 Q(K,k') 
2tt 2 

b 



E{k)-^f+iQ 



(22) 



Here, the scattering energy E(k) is taken relative to the 
energy at the continuum edge E cont (K). The in- medium 
scattering phase shift follows according to eq.(^Tj) and 

eq.© 

1 -Re J (K;E cont (K)+E(k)) 



cot S(K;E(k)) = 



Im J(K;E coat {K)+E(k)) 



(23) 



The small k expansion (g) can be applied to the in- 
medium phase shifts S in the same way as for the free 
phase shift and provides now a quantity which is a direct 
generalization of the free scattering length ao- In anal- 
ogy to the phase shift it is named in-medium scattering 
length a(/i, T) and follows in the center-of-mass system 
(K = 0) from 



a(/x, T) = — lim 



ImJ(^,^ cont (0) + g(fc)) 
k"'o k [1 - Re J(K, E cont (0) + E(k))} 



(24) 



If we consider a repulsive potential (no bound states) and 
cooled spin-polarized atoms (triplet scattering only) as 
used in magnetic traps, the generalized Beth-Uhlenbeck 
formula eq. ( |l4| ) simplifies after partial integration to 



" 1 ./ii 

8{K-E) 



dKK 2 
2tt 2 



dE 



2kBT[cosh((E + E cont {K) - 2^*)/k B T) - 1] 



(25) 



Now we focus on the case of a contact potential. This 
is done to compare with classical results as eq.(Q) which 
are also obtained with a zero range pseudopotential. A 
physical reason is the assumption <C ao ^ A which 
makes the first term in expansion (|^) the dominant one. 
It is valid for the atomic systems considered below if we 
have low temperatures and a potential range r in the 
order of the atomic hard core. 

For a zero range potential (/3 — > oo in cq.([i6|)) the 
T-matrix takes the form 



( 



Ti 2 k 2 \ AttH^ a 



T K:E = 



X (l + ia k+^ I dy^L{Q{K,Jy)-l)) \ (26) 



y 



Here, the free T-matrix is first integrated out and then 
the limit f3 — > oo is taken p^j. This procedure avoids 
the introduction of a cut-off parameter to regularize the 
integration. 
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The in-medium scattering phase shift for the zero range 
interaction reads 



S { K:E = 



arctan 



h 2 k 2 



aok Q{K, k) 



-l + ffdy-^^y(Q(K,^y)-l), 



(27) 



It differs from the free scattering phase, which is sim- 
ply So(k) = — arctan (do k), by the occurrence of the Bose 
enhancement factor Q according to eq. ( pl| ) in the nom- 
inator and in the principal value integration in the de- 
nominator. Then, the in-medium scattering length for 
the zero range interaction is 



t(/i, T) = a \ 



l + 29i(0) 



.1 + ^ j dk' 9l {e k ,), 




(28) 



The free scattering length ao is modified by the Bose 
functions in the nominator and denominator. For ji* — > 
the Bose pole in the nominator dominates the expression 
( p8| ) and the in-medium scattering length diverges as one 
approaches the onset of Bosc-Einstein condensation. 

To apply the considerations above we focus on two sys- 
tems which play an important role in the current discus- 
sion on atomic Bose gases. For the first one, sodium 23 Na, 
a free ground state scattering length of ±(92 ± 25)as (in 
units of the Bohr radius as) is obtained from cross sec- 
tion measurement Spectroscopic data are in agree- 
ment with 64 < ao/a^ < 152 [|l2|, i.e. repulsive interac- 
tion is favoured and bound state formation is excluded. 
For the second example, rubidium Rb, the ground state 
triplet scattering length ao ~ 5 nm is of the same order 
as for sodium M. Spectra are analyzed to yield values 
85 < a /a B < 140 @. For our exploratory calculation 
we have chosen a free scattering length a Q = 92 as for 
sodium and ao = 95 as for rubidium 87 Rb to parametrize 
the T- matrix (Eq). 



V. RESULTS AND DISCUSSION 

Figure [l] shows the result for the in-medium scatter- 
ing phase shift of sodium according to eq.(27). For 
some densities the dependence on relative momentum 
is plotted and compared to the free scattering phase 
shift. The temperature T = 5/iK has been chosen to 
have similar values as reached in the recent experiments. 
E.g., in the sodium experiment densities exceeding 



10 cm" J have been produced and a condensate has 
been observed at temperatures below 2 /iK. The den- 
sities in our calculation are below the critical density of 
^bec — 6 x 10 cm" 3 at T — 5/iK. In general, medium 



modifications of the phase shift compared to the free 
scattering phase shift become small (< 1%) for scatter- 
ing energies h 2 k 2 /m > (10...20)/cbT. The influence of 
the medium enlarges the absolut value of the phase shift 
and leads to an extra minimum if one increases the den- 
sity towards the critical value. Furthermore, we see that 
the medium increases the derivative of the phase shift 
at small momenta. This will result in an increase of the 
scattering length as demonstrated in figure ||. Since the 
free scattering length appears to be of the same size for 
sodium and rubidium the major difference in the results 
originates from the different atomic mass. Consequently, 
compared to sodium at the same temperature a similar 
behaviour of scattering shift and length occurs for rubid- 
ium at densities about seven to eight times higher. 
Figure || shows that the scattering length according to 
eq.(|2|) enlarges with increasing density due to the Bose 
enhancement. We have plotted the density dependence 
of the in-medium scattering length of sodium for two dif- 
ferent temperatures. It shows a strong increase if the 
density approaches the critical density tibec- Accord- 
ing to our model calculation for sodium the value of the 
scattering length is doubled for T = 5/iK at a density of 
1 x 10 14 cm- 3 (n BE c = 6x 10 14 cm~ 3 ) and for T = 50/iK 
at a density of 3.7 x 10 15 cm~ 3 (n B EC = 1.9 x 10 16 cm" 3 ). 
For rubidium 87 Rb at T = 5/iK the scattering length is 
doubled at 8 x 10 14 cm- 3 (n BE c = 4.5 x 10 15 cm~ 3 ). If 
one measures the cross section a clear enlargement of the 
scattering length should be observed already at densities 
about 5 times smaller than the critical density for Bose- 
Einstein condensation. 

In figure || we have evaluated the generalized Beth- 
Uhlenbeck equation as given in eq. (p5|) and have com- 
pared the non-ideality coefficient B*(/j,, T) of eq.([l5|) with 
the classical second virial coefficient B n {T) of eq.(Q). 
The relation of both is plotted for 23 Na and 87 Rb ver- 
sus the free density of quasi-bosons at a fixed temper- 
ature T = 5/iK. Approaching the critical density, the 
non-ideality coefficient B*(fi,T) increases strongly. The 
mass difference causes the respective behaviour for rubid- 
ium to occur at roughly seven times higher densities then 
for sodium at this temperature. In the low density limit 
B* approaches the second virial coefficient B n . Figure 
^ shows a minor deviation especially for 87 Rb from the 
expected low density limit B* /B n = 1. It can be traced 
back to the limited validity of the condition ao <C A nec- 
essary to derive (||) . The classical second virial coefficient 
B n as the low-density limit is slightly overestimated for 
higher temperatures and greater masses. We stress that 
the absolute value of the non-ideality coefficient B* (/i, T) 
is monotonically increasing with density. 

We have also performed calculations with a separable 
potential of finite range. The model parameter Vq and 
(3 of the potential ([l6]) have been fitted to the free scat- 
tering length ao and the potential range Vq. The latter 
one characterizes the short range behaviour of the inter- 
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action. Since no exact datas are available, we have taken 
a value in the order of the hard core radius of the atoms. 
A theoretical value of the hard core radius is about Aclb 
(from [^3| at higher temperatures). The results with this 
finite range interaction have been found to be not sen- 
sitive if compared to the results presented in the figures 
1-3 for zero range interaction. This is not very surprising 
since ro is much smaller than the free scattering length 
do. But we are cautious to state that the results are 
in general independend of the potential range since we 
know neither ro nor a hard core radius at this very low 
temperatures. 

We remind that we restricted our calculations to the 
normal phase. Improving the ordinary quasiparticle 
picture, correlations are included within the so-called 
extended quasiparticle approximation. In the Bosc- 
condensed phase the common approaches within Hartree- 
Fock-Bogoliubov approximation (see HHI) should be 
improved on the same level by including correlations be- 
yond mcanficld. 

Furthermore, our calculations are restricted to the case 
of repulsive interacting bosons. The consideration of an 
attractive interaction leads to the possibility of an ad- 
ditional liquid-gas phase transition and possibly to the 
formation of bound states. The question of a superposi- 
tion of a gas-liquid or gas-solid phase transition and the 
Bose-Einstein condensation in an attractive atomic Bose 
gas like lithium 7 Li is discussed in j2j|. There, the equa- 
tion of state is extended to the Bose-condensed region. In 
the normal phase it is considered in the approximation of 
eq.(|^) which does not include in- medium scattering and 
possible contributions from bound states. A treatment in 
the framework of the generalized Beth-Uhlenbeck equa- 
tion ( pT| ) would provide a better understanding of the 
possible phase transitions. 

In a generalization of the derived results to non- 
equilibrium processes, transport coefficients should be 
modified not only due to final state Bose enhancement 
, but as well due to intrinsic medium effects. The lad- 
der T-matrix should serve as an ingredient for more ad- 
vanced formulations of kinetic equations including strong 
correlations H]. 

In conclusion we suggest that in the recent experiments 
with Bose gases, such as sodium and rubidium in atom 
traps, density and temperature regions are within reach 
where the scattering properties and the thermodynamic 
properties are influenced by the medium. We expect that 
the medium modifications have direct impact on the ex- 
perimentally extracted values of the cross section and 
hence the scattering length. 
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FIG. 1. Free and in-medium scattering phase shift versus 
decadic logarithm of relative momentum for sodium. Tem- 
perature is T = 5^iK. 



FIG. 3. The nonideality coefficient B* (//, T) in units of the 
classical second virial coefficient B n (T) for sodium and rubid- 
ium versus free density rif ree (/u, T) at temperature T = 5/iK. 
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FIG. 2. The in-medium scattering length a(fi, T) in units of 
the free scattering length do versus free density for sodium at 
temperatures T = 5/iK and T = 50/itK. The critical densities 
are marked. 
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